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Abstract 

Recent results on top quark physics with the D0 experiment in pp collisions at 
\fs = 1.8 TeV for an integrated luminosity of 125 pb _1 are reported. The direct 
measurement of the top quark mass uses single lepton and dilepton events, giving 
the result rritop = 172.1 ± 7.1 GeV/c 2 . The measurement of the tt production cross 
section includes analyses from 9 top decay channels: dilepton (tt — ► e/i, ee, and /i/i), 
electron and neutrino (tt — > ev), single leptons (tt — > e + jets, tt — > la + jets) with 
and without b tagging, and all-jets (tt — * 6 jets). We measure the tt production 
cross section to be 5.9± 1.7 pb at rritop = 172.1 GeV/c 2 . Combining the D0 and the 
CDF measurements of the top quark mass and combining their tt cross sections, in 
both cases taking into account error correlations, yields unofficially a t j = 6.7 ± 1.3 
pb at an averaged top quark mass of rritop = 173.8±5.2 GeV/c 2 . Preliminary results 
on a search for charged Higgs production in top events, t — ► H + b, are presented. 



presented at the 1998 Rencontres de Physique de la Vallee d'Aoste on Results and Perspectives in 
Particle Physics, La Thuile, Italy, March 1-7, 1998. 



1 Introduction 



Since the top quark was discovered |T], |2[ by the D0 and CDF experiments at the Fermilab 
Tevatron collider in 1995, much of the effort has been focused on using all available 
information on the cadidate events to measure the it production cross section and the 
top quark mass as precisely as possible. 

In the Tevatron pp collider at a center of mass energy of 1.8 TeV top quarks are 
predominantly produced in pairs through qq annihilation (~ 90%) or gluon fusion (~ 
10%). According to the standard model a top quark decays almost 100% of the time into 
a W and a b quark. Each W boson decays either into a charged lepton and a neutrino 
or into a pair of quarks. Our analysis channels are classified on the basis of the W boson 
decay. The cleanest channel is the dilepton channel, it — > 2 leptons + 2 neutrinos + 2 
jets, where both W bosons decay into leptons. The branching ratio is ^ for it decays into 
e/i, ee, and \i\x channels. The single lepton (e or \i) decay channel, where one W boson 
decays leptonically and the other W boson decays hadronically, has a branching ratio of 
|y, but it has a sizable background from W^+jets production. The all-jets decay channel, 
where both W bosons decay hadronically, has a branching ratio of ||, but it has a huge 
QCD multijet background. 

This paper reports the recent D0 results on the it production cross section using the 
dilepton, single lepton [[J, and all-jets 0] decay channels, and on the direct measurement 
of the top quark mass using the single lepton 0, |J and dilepton events [0]. 

The results reported in this paper use the entire data sample with an integrated 
luminosity of about 125 pb _1 , collected during the 1992-1996 collider run. Since our 
report on the discovery of the top quark fjl|], our data sample has more than doubled. D0 
has optimized the analysis to maximize the expected precision of the it production cross 
section measurement and improved the techniques on the measurements of the top quark 
mass. 

2 Measurement of it Production Cross Section 

The triggers and the particle identification algorithms used in these analyses are described 
in detail in D0 publications [p], §]. 
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2.1 Dilepton Decay Channels 

The signature of it events in the dilepton decay channels is two high p T isolated leptons, 
two or more jets and large missing ($t ) due to the presence of the two neutrinos. 

The offline event selection requires that the isolated electrons have E? > 20 GeV with 
I77I < 2.5 and that the isolated muons have pt > 15 GeV with \r]\ < 1.7. At least two jets 
are required to be reconstructed with a transverse energy above 20 GeV with \rj\ < 2.5. 
All jets are reconstructed using a cone algorithm with radius 0.5 in rj-4> space. The ]£ T is 
required to be above 25 GeV for the ee channel and to be above 20 GeV for the e\i 
channel. In addition, we apply a cut on a variable Ht, which is defined as the scalar sum 
of the Ex's of all the jets in the event plus the leading electron Ef. We expect a higher 
Ht for it events than for background events. We require the Ht to be above 120 GeV 
for the ee and e/x channels, and to be above 100 GeV for the \i\x channel. 

The background in these channels is mainly from Z° decays, Drell-Yan, vector boson 
pair events and events with misidentified leptons. 

Five events passed the above selection criteria for the dilepton decay channel, with 
estimated background 1.4 ± 0.4 events (see Table [I]). 

2.2 Single Lepton Decay Channels 

The signature of it events in the single lepton channels is a high pt isolated lepton, 
large $ T an d high jet multiplicity. We divide the single lepton channel analysis into 
two complementary and orthogonal analyses: an event shape analysis and a b tagging 
analysis. The main sources of background in the single lepton channels are W + jet 
events with high jet multiplicity and multijet QCD events with fake leptons and $ T due 
to mismeasurement. 

In the event shape analysis the it — > e + jets and t i — > \i + jets events are selected 
by use of topological and kinematic cuts. The offline event selection requires one isolated 
lepton with E T > 20 GeV with \rj e \ < 2.0 or |^ t | < 1.7 and at least 4 jets with E T > 
15 GeV with \r)jet\ < 2.0. The $ T requirement is 25 GeV for the e + jets channel 
and 20 GeV for the \x + jets channel. The event shape analysis uses the global event 
variable aplanarity, A, which is defined as | of the smallest normalized eigenvalue of the 3- 
momentum tensor of the W boson and jet momenta in the laboratory frame. Background 
events typically have smaller A than signal events. Another variable used in the event 
shape analysis is Ht, which is defined as the sum of the Et's of all the jets in the event. 
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The it events are expected to have a higher Ht than W^+jets events. In addition we also 
require the total leptonic E T , E^, which is defined as the sum of the charged lepton E T 
and $ T , to be above 60 GeV to reject the multijet QCD background. We require the 
pseudorapidity r]w of the W boson that decays leptonically to have \rjw\ < 2.0 to obtain 
better agreement between background control samples from data and the W+jets Monte 
Carlo (MC) samples §. 
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Figure 1: Aplanarity vs Ht for single lepton data, it MC events, multijet background and 
W+4 jets MC background. The dashed lines indicate the cuts. 

Figure [i] shows scatter plots of A vs H T for D0 single lepton data, it (m t o P = 
170 GeV/c 2 ) MC events, and the two background sources: multijet background and W + 
4 jets VECBOS MC events. Based on our optimization procedure using it MC events, 
we define the region A > 0.065 and Ht > 180 GeV as the signal region. 

W + jets background events and the multijet QCD events are the main source of 
backgrounds in the lepton + jets channel. Since the number of W + jets events is 
expected to decrease exponentially as a function of the jet multiplicity [fTofl , by fitting the 
number of W + jets events at the lower jet multiplicity and extrapolating it to high jet 
multiplicities, we can estimate the number of W + jets background events in the data 
sample. We estimate the QCD background from the data itself using the measured jet 
misidentification probability. 

Nineteen events survived all the cuts in the event shape analysis for the single lepton 
channel, with estimated background 8.7 ± 1.7 events (see Table [1|). 

In the b tagging analysis, the background for the single lepton channel can be sig- 
nificantly reduced by requiring that one of the jets is tagged as a 6-jet. We tag 6's by 
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detecting a muon in a jet. About 20% of it events have a detectable \i in a jet compared 
to only about 2% of the W + (>)3jets background events. A tag muon is required to have 
p T > 4 GeV and the distance AR between the muon and a jet in the (p-rj plane must be 
less than 0.5. 

The offline event selection in the b tagging analysis requires one isolated lepton with 
E T > 20 GeV with \rj e \ < 2.0 or |r^| < 1.7, $ T > 20 GeV, and 3 or more jets with 
Et > 20 GeV with \i]j et \ < 2.0. Since we require a tag muon in the event, we use 
looser cuts on aplanarity and Ht compared to the event shape analysis: A > 0.040 and 
H T > 110 GeV. Figure § shows the distribution of the jet multiplicity for background 
events and for single lepton data before the A and H T cuts. The data agrees with the 
number of events from W + jets and QCD processes in the low jet multiplicity region. 
For the high multiplicity bin we can see a clear excess above background even without 
the A and H T cuts. 
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Figure 2: Jet multiplicity distribution for single lepton fc-tag events (before applying A 
and Ht cuts), compared to background estimates. 

The main background in the b tagging analysis is also M^+jets production. The back- 
ground is estimated by applying the muon tag rate, which is determined by the fraction 
of jets tagged in multijet events, to the jets in the background sample after all other cuts 
described above are applied. 

Eleven events passed the above selection criteria in the b tagging analysis of the single 
lepton channel, with estimated background 2.4 ±0.5 events (see Table |l|). 

2.3 ev Decay Channel 

In order to identify it events in which neutrinos from both W decays carry much of the 
transverse momentum, we perform an ev decay channel analysis. This analysis is mainly 
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focused on selecting those tt events that fail the event selection criteria for the ee, e\i and 
e+jets decay channels. The signature of the tt events in the ev decay channel is one high 
Pt isolated electron, two or more jets and very large lfi T which together with the electron, 
forms a transverse mass much higher than the W mass. 

The event selection requires events with an isolated electron with Et > 20 GeV with 
|?7 e | < 1.1 and at least two jets with Et > 30 GeV with \r)j et \ < 2.0. We require that the 
]p T is above 50 GeV and the transverse mass of the $ T and the electron is above 115 GeV. 
The dominant background in this channel is VT+jets production, W pair production and 
a misidentified electron with $ T in QCD multijet events. 

Four events survived the selection criteria for the ev decay channel, with estimated 
bachground 1.2 ± 0.4 events (see Table D). 

2.4 All-jets Decay Channel 

The decay of the top quark through the multi-jets channel represents the dominant decay 
mode in the standard model (44%). But, without the signature of an isolated lepton 
in the event, the QCD multi-jets background overwhelms the signal by roughly 1:2000 
before any other cuts are applied. This has made the search for top in this channel very 
challenging. Significant gains in this channel have been made by employing multivariate 
techniques and a number of novel techniques which are optimized for this search. The 
current analysis [|J includes the full Run 1 data set. 

Artificial neural networks constitute a powerful nonlinear extension of conventional 
methods of mult i- dimensional data analysis ||11|| . They have been proven extremely well 
suited to the multi-jets search because they handle information from a large number of 
inputs and properly account for the nonlinear correlations between these inputs. We note 
that the output of the neural net is simply a mapping between the multi-dimensional space 
of variables describing our events and a one- dimensional output space. The analysis may 
proceed by setting a threshold on the neural net output, corresponding to a cut surface 
in multi-dimensional space. 

This analysis relies on a set of simple topological variables which discriminate between 
the fundamental properties of the tt signal and the QCD background. For the background, 
these properties are: 

• The events have a lower overall energy scale; leading jets have lower Et', multi-jet 
invariant masses are small. 
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• The additional radiated jets are soft (have low Et)- 

• The event topology is more planar (less spherical than top). 

• The jets are more forward-backward in rapidity (less central than top). 

We employ 10 kinematic parameters which embody these properties. In addition to 
the kinematic parameters, we use three other parameters based on jet widths, px of the 
tagging muon and the mass likelihood. This set of 13 variables was used as the set of 



inputs to the neural net. We used the back-propagation learning algorithm in Jetnet [12] 
to train the network and provide a single variable NN (the neural network output in the 
range (0,1)). 

At the simplest level, each t(i) decays into a b(b) and W + (W~) boson, and the W 
bosons into light quarks. Barring extra gluon bremsstrahlung, this represents six quark- 
jets in the final state. The average jet multiplicity for Herwig ti events (mt=180 GeV/c 2 ) 
is 6.9, implying that the contribution from gluons is relatively small. Conversely, the QCD 
multi-jet background has final-state jets that originate predominantly from underlying 
gluon radiation. Although gluon splitting can take place, producing both quark and 
gluon jets, it is expected that there will be a large enhancement in gluon content for QCD 
jet production. A Fisher discriminant, based on the difference between quark and gluon 
jet widths is used to differentiate between ti signal and QCD background. 

The pt of the tagging muon helps to further differentiate between signal and QCD 
background. Not only does the fragmentation of b quarks produce higher px objects, but 
the b itself is more boosted in top events than in background. Thus, the muon px spectrum 
is significantly harder in ti events. The muon px serves as a useful tool in differentiating 
between signal and background. 

A mass log-likelihood variable provides discrimination between signal and background 
by requiring that there are combinations of jets consistent with being from two top decays 
and the subsequent W — > qq. The presence of two W's in it events provides significant 
rejection against QCD background. A further requirement that the two reconstructed top 
quarks have equivalent masses provides additional discrimination. 

The background NN distribution is based upon data events, selected by identical 
criteria as the signal sample except that a tagging muon is not required. To these events, 
a muon is added with probability and kinematic distributions as observed in our data. 
The it NN distributions are taken from the Herwig MC, and checked with Isajet. 
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Figure 3: Fit to the distributions in neural net output for signal and background. The normal- 
izations of each are fitted to the data (diamonds) . 



Once the final neural net has been trained, the relative amount of signal and back- 
ground in the data may be determined by a fit for the amounts for signal and background. 
Because the shapes of the expected signal and background distributions in neural net out- 
put differ significantly, the relative amounts of each may be disentangled. This serves as 
a useful check of the overall tag rate function and simultaneously provides the top cross 
section. This technique benefits from the fact that additional information on the top cross 
section is extracted even in regions where the background level is large. The data (dia- 
monds), and the expected distribution for background and for signal are shown in Fig. 
along with the combined fit. The result is, <j t t = 7.1 ± 2.8 ± 1.5 pb at m t = 172.1 GeV/c 2 , 
in good agreement with the D0 result for the leptonic channels ||. Applying a cut on the 
neural net output of 0.85 (events shown in Table [[]), and using the absolute background 
predicted, yields a consistent cross section. Additionally, a more restrictive cut on the 
neural net output of 0.94, chosen for maximum expected significance using Monte Carlo 
data, leaves 18 events with an expected background of 6.9 ± 0.9. The probability of an 
upward fluctuation of the background to the observed signal is equivalent to 3.2 standard 
deviations. 



2.5 tt Production Cross Section 

Table [1] shows the number of observed events for all nine decay channels, the esti- 
mated number of background events, and the expected number of tt events for m top = 
170 GeV/c 2 (tt production cross sections using ref. |P~3"[). 

The top production cross section in the leptonic channels is calculated using the for- 

mula o t i = — 1 — - where N{ is the number of observed events, Bi is the expected 
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Table 1: Event yields 



channel 


events 


estimated 


expected signal 




observed 


background 


m t =170 GeV/c 2 


ee 


1 


0.5±0.1 


1.2±0.2 


e// 


3 


0.2±0.2 


2.2±0.5 


nn 


1 


0.7±0.3 


0.6±0.1 


e+jets(event shape) 


9 


4.5±0.9 


8.6±1.6 


^±jets(event shape) 


10 


4.2±1.0 


5.5±1.7 


e+jets(6-tag) 


5 


1.1±0.4 


3.6±0.6 


/i+jets(6-tag) 


6 


1.4±0.2 


2.3±0.5 


eu 


4 


1.2±0.4 


1.7±0.5 


total (leptonic channels) 


39 


13.7±2.2 


25.7±4.6 


all-jets 


41 


24.8±2.4 


15.9±2.6 



background, £j is the detection efficiency for top, Bi is the branching ratio and Lj is in- 
tegrated luminosity for channel i, where i=l,...,8. This cross section is then combined 
with the aforementioned cross section measured in the all-jets decay channel taking into 
account appropriate correlations. 
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Figure 4: Top cross section vs top quark mass 

Figure | shows the D0 tt cross section as a function of top quark mass compared with 
theoretical predictions. We measure the it production cross section to be a t j = 5.9 ±1.7 
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pb at our measured top quark mass of m tov = 172.1 GeV/c 2 (see section |) in good 
agreement with the standard model predictions [T3|, [15|, [16|] . 
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Figure 5: Summary of tt cross section measurements and theoretical calculations 



Combining the D0 and the CDF |L7j tt measured cross sections taking into account 
correlations between systematic errors yields unofficially a t j = 6.7 ± 1.3 pb at an averaged 
top quark mass of m top = 174 GeV/c 2 . Figure [5] shows the D0 and the CDF measured 
tt cross sections as well as the unofficial Tevatron average and a theoretical band which 
contains the four calculations. 



3 Measurement of the Top Quark Mass 

Direct measurements of the top quark mass using single lepton and dilepton events are 
reported here. 



3.1 Single Lepton Events 
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3.1.1 Event Selection 



The initial event selection for the top quark mass measurement analysis is similar to that 
used in the cross section analysis for the single lepton channel, but without A and 
cuts. We require one isolated lepton, e or /x, with E T > 20 GeV and \r) e \ < 2.0 or 
1^1 < 1.7. We also require $ T > 25 GeV for e + jets and 20 GeV for ji + jets. Only the 
events with four or more jets having Et > 15 GeV with |?7j et | < 2.0 are used. To suppress 
background in untagged events, we require E% > 60 GeV and \r)w\ < 2 for the W — > Iv. 

Ninety events passed the event selection requirements. Among these, seven events are 
6-tagged events. 

3.1.2 Fitting Algorithm 

For each event passing the above selection cuts, we make a two constraint kinematic 
fit |TJ|] to the it — > £ + jets hypothesis by minimizing a x 2 = ( v ~ v*) T G(v — v*), where 
v(v*) is the vector of the measured (fit) variables and G _1 is its error matrix ||. Both 
reconstructed W masses are constrained to equal the W mass and we assume both t and 
t quarks have the same fit mass, m^t- Kinematic fits were performed on all permutations 
of the jet assignments of the four highest Et jets, with the provision that muon-tagged 
jets were always assigned to a 6-quark in the fit. Each fit yields a fitted mass value, mg t 
and a x 2 - The fit with the lowest \ 2 is chosen to describe the event. Only the events with 
X 2 < 10 are used in the top quark mass determination. 77 events passed the x 2 cu ^ and 
among them five events are 6-tagged events. Although mg t is strongly correlated with 
the top quark mass, m top , m fit is not the same as m top , because of gluon radiation and 
permutation ambiguities. 

3.1.3 Mass and Error Determination 

To further separate the signal and background events we use variables that can provide 
good separation between it and background events without much correlation to the fit- 
ted mass. The following four variables are chosen to compute the top quark likelihood 
discriminant (D): 

• $ T . 

• A. 

• H T 2/^\pz\, where H T 2 is the H T without the E T of the leading jet. 
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minimum AR between all pairs of the 
jets and E™ m is the smaller jet from the minimum AR pair. 



D0 uses two discriminants ||. One is the Dlb (low bias) method, in which we parametrize 
Ci{xi) = Si(xi)/bi(xi), where s, and bi are the top signal and background densities in each 
variable, integrating over the others. We then form the log likelihood ln£ = J2i &i hi A, 
where the weights U{ are adjusted slightly away from unity to nullify the average correla- 
tion ("bias") of C with mg t . For each event we set D LB = £/(l + C). The data are then 
divided into two bins: a low signal-to-noise bin and a high signal-to-noise bin, according 
to whether the LB cut is passed. The LB cut is passed if either Dlb > 0.43 and Hti > 90 



GeV, or if a b tag exists in the event. Another method uses a neural network [19| with 
the same four variables as input, five hidden nodes, and one node with output D^n- We 
divide the data into ten bins in D NN . Figure |] shows the distribution of the discriminants 
-D L b an d -D NN for signal and background. They indicate that either discriminant provides 
good discrimination. 
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Figure 6: Distribution of discriminants, D LB and D NN for 175 GeV/c 2 top MC (shaded 
histogram) and for background (unshaded histogram). 

Since the selected event sample contains both it and background events, we make a 
two dimensional likelihood fit of the event sample to the sum of the expected it signal 
plus background in the mf it vs. D plane. We make an independent likelihood fit for 
each top quark mass hypothesis. We use the herwig |20| MC to simulate it events. We 
estimate background using a combination of P^+jets events from the vecbos MC and 
fake lepton events obtained directly from D0 data. 
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Figures [7] (a) and (b) show the distributions of mat for data (a) passing and (b) failing 
the LB cut. The histograms are data, the dots are the predicted mixture of signal plus 
background, and triangles are background. Figure (c) shows the log of the fit likelihood 
L vs. true top quark mass mt op for both the LB and NN fits. The curves are quadratic 
fits to the lowest point and its 8 nearest neighbor points. The minimum position of 
each curve yields the measured top quark mass. The width of the curve at 0.5 above 
the minimum determines the statistical error of the measurement. The LB fit yields 
m top = 174.0 ± 5.6 (stat.) GeV/c 2 . The NN fit yields m top = 171.3 ± 6.0 (stat.) GeV/c 2 . 

There are several sources of systematic error |5|, ^j. The major uncertainties are from 
the jet energy scale and the MC modelling. We assign a jet energy scale error of ±(2.5% + 
0.5 GeV) based on a detailed study of 7+jet events in data and MC, particularly focused 
on the dependence of the Et balance upon rj of the jet, and checked by the Et balance 
in Z±jet events. This leads to an error on m top of ±4.0 GeV/c 2 . The uncertainties in 
the MC modeling of top and W+jet production are estimated to have a systematic error 
on the top mass of ±3.1 GeV/c 2 . Other effects, including calorimeter noise, multiple pp 
interactions, and differences in fits to InL, contribute ±2.2 GeV/c 2 . All systematic errors 
sum in quadrature to ±5.5 GeV/c 2 . 
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Figure 7: (a-b) Fit top quark mass distribution for events (a) passing or (b) failing the 
LB cut. (c) Log of likelihood L vs. true top quark mass mt op for both LB and NN fits. 



Combining m top from both methods, LB and NN, we determine the top quark mass to 
be m top = 173.3 ± 5.6 (stat) ± 5.5 (syst) GeV/c 2 or m top = 173.3 ± 7.8 GeV/c 2 , allowing 
for the (88 ±4)% correlation between two methods. 
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Figure 8: Sum of the normalized candidate weights grouped into the five bins considered 
in the maximum likelihood fit (circles) for the D0 MWT (left) and i/WT (right) analyses. 
The uncertainty on these points is taken from the RMS spread of the weights in MC 
studies. Also shown are the average weights from the best-fit background (dashed) and 
signal plus background (solid). The -In L distributions and quadratic fits are inset. 



3.2 Dilepton Events 

Due to the presence of two neutrinos, dilepton events do not contain enough information 
for a constrained fit. Therefore, we resort to using mass estimators other than the recon- 
structed mass. The same technique that was used for the lepton plus jets mass analyses is 
used to calculate a top quark mass likelihood in the dilepton mass analysis. D0 supplies 
the missing constraint in the problem by assuming a top quark mass and reconstructing 
the event for every assumed mass 0. Based on the reconstructed final state a weight is 
computed which characterizes the probablity for this event to be from a it decay with 
the assumed mass. There are two algorithms used to determine the weight. The matrix 
element weighting (MWT) method uses the proton structure functions and the probabil- 
ity density function for the energy of the decay lepton in the rest frame of the top quark 
(an extension of Ref. ||22|| ). The neutrino weighting method (z/WT) assigns the weight 
based on the available phase space for the neutrinos, consistent with the measured lfi T . A 
maximum likelihood fit is performed to the shape of the normalized weight curve summed 
over all 6 dilepton mass events using MC derived probability density functions for signal 
and background (see Fig. |8|). The results for the two analyses are in excellent agreement 
: m t = 168.1 ± 12.4(stat) GeV/c 2 (MWT) and m t = 169.9 ± 14.8(stat) GeV/c 2 (z/WT) 
with a systematic error of 3.7 GeV/c 2 . By combining the two results, taking into account 
the correlations (77%), we determine the top quark mass in the dilepton channels to be 
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m t = 168.4 ± 12.3(stat) ± 3.6(sys) GeV/c 2 or m t = 168.4 ± 12.8 GeV/c 2 . 

D0 has also combined its dilepton and lepton+jets mass measurements into a single 
top quark mass measurement of m t = 172.1 ± 5.2(stat) ± 4.9(sys) GeV/c 2 or m t = 
172.1 ±7.1 GeV/c 2 . 



3.3 Combined Top Quark Mass 



CDF dilepton (1) 
CDF lepton+jets (2) 
CDF all jets (3) 

DO dilepton (4) 
DO lepton+jets (5) 
(DO combined) 

Tevatron lepton+jets 
Tevatron Average 




161.0 ± 20.0 
175.9 ± 6.9 
186.0 ± 15.2 
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174.8 ± 5.5 
173.8 ± 5.2 
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Figure 9: Summary of direct measurements of the top quark mass 



Combining the D0 and the CDF |T7j top quark mass measurements taking into ac- 
count correlations between systematic errors yield unofficially mt op = 173.8 ±5.2 GeV/c 2 . 
Figure |9] shows the D0 and the CDF measured top quark mass as well as the unofficial 
Tevatron averages. 



4 Charged Higgs Search 
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4.1 Introduction and Strategy 



The "minimal" Higgs sector of the Standard Model (SM) is comprised of only one complex 
Higgs doublet, which manifests itself in a single physical neutral Higgs scalar H° (m#o is 
a free parameter of the SM). However, the scalar sector of the SM is not yet constrained 
by experiment, and it is possible that the Higgs content is richer and its structure more 
complicated. A strong motivation for extending the Higgs sector comes from supersym- 
metry, where two Higgs doublets are needed to provide masses to the charged fermions. 
Two-higgs doublets appear in many theories, both in the context of the SM and extended 
models such as the E 6 superstring-inspired model. Charged Higgs bosons appear in these 
as well as in several more general theories [|23 . 



The minimal extension to the Higgs sector has two complex SU(2)l (isospin T3) Higgs 
doublets with opposite hypercharge [23 1. In the Minimal Supersymmetric Model (MSSM) 
there are three parameters: the mass of the W boson (raw), the ratio of the two neutral 
Higgs VEVs (tan ft), and the mass of the charged scalar (m H ±). If kinematically allowed, 
BR(t — > H + b ) can be comparable to, or even larger than, BR(t — > W + b ) depending on 
m^+and tan ft. At the Tevatron one can therefore search for a light charged Higgs boson 
in the decay of top quarks. 

There is a direct lower limit from LEP experiments corresponding to m#± > 75 GeV 
at 95% CL p5|. Based on their direct and indirect searches, the CDF experiment has 



recently excluded additional regions in the [mjj+ , tan ft] plane P5[ . Also, an independent 
analysis of CDF data led to similar conclusions for high tan/3[^7|. Our analysis is based 
on a somewhat different strategy and is sensitive to similar but not identical regions of 
the parameter space. 

The value of the it production cross section and the mass of the top quark from direct 
observation at D0 || || are based on the assumption of BR(t — > W + b) = 1. The topic 
of present interest is based on the violation of this assumption. We therefore cannot use 
those measurements in this analysis. Consequently, in the process of setting limits in the 
[m,H+, tan/?] plane, a (it) and m t serve as input parameters. Production of it takes place 
via strong interactions, and the cross section is not affected by the existence of a charged 
Higgs boson lighter than the top quark. One can therefore use a theoretical prediction 
based purely on the SM for a(ti). There is also no rigorous reason to use the directly 
measured value of m t when considering a region where BR(t — > H + b) « 1, but in the 
absence of a compelling reason to start elsewhere, one can choose m t in the 170 - 175 
GeV range as a starting point. In our analysis we have considered so far charged Higgs 
predominantly decaying into the heaviest quarks or leptons, i.e. H + — > cs or H + — > t + u t , 
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(and the charged conjugate decays). The branching ratios for these decays depend on 
taxi p. Neglecting other decays into e + h , e ,fi + u^ and ud, BR(H — > cs)+BR(H — > tu t ) = \. 

The selection criteria of the D0 standard measurement of a(ti) were optimized for 
ti — > W + bW~b. In the end, the number of observed events in excess of the estimated 
backgrounds is interpreted as arising entirely from ti — > W + bW~b, and a{ti) was calcu- 
lated assuming BR(t — > W + b) = 1. In our indirect search, we replace this last assumption 
with a more general one: BR(t — > W + b) + BR(t — > H + b) = 1. We then appply the same 
selection criteria to our simulated signal and determine, for given values of , m t and 
<j{tt), what values of tan/3 are unlikely to produce the number of events observed in our 
data sample (n obs ). With selection criteria tuned for high acceptance to ti — > W + bW~b 
events, the value of n b. s — (nbkg) is statistically significant, and in excellent agreement with 
the predicted a (it) for BR(t — > W + b) « 1 (i.e., BR(t — > H + b) « 0). Therefore, regions 
of parameter space that imply a high BR(t — > H + b) can be excluded if the acceptance 
(using our normal selection criteria) for ti — * H + X is small, because it is not possible 
to explain the large value of n obs — (ribkg) in these regions. This is the strategy of our 
"disappearance" analysis. 

For the time being, we limit ourselves to the lepton+jets final states, which accounts for 
30 of the 39 ti candidate events on which the current D0 standard measurement of a(ti) 
is based. One very appealing feature of the indirect search is that the selection criteria are 
already tuned and the non-top backgrounds, as well as the systematic uncertainties, have 
been thoroughly studied and are well understood through the previous analysis. For an 
indirect search, one generates the charged Higgs signal using Monte Carlo, and studies its 
acceptance for the standard event selection criteria. Events are generated for one channel 
at a time for any given m#+. Subsequently, a sample for a particular value of tan/5can 
be obtained by adding events from different channels in a proportion appropriate for that 
[m H +, tan/5]. A single set of selection criteria is applied to all samples, each of which 
correspond to a unique [m H +, tan/?]. 

4.2 Results and Discussion 

For a given set of the three parameters (m t , m#+, a(ti)), we determine the probability 
distribution of tan/3. There is no generally applicable strict upper or lower limit on tan/3. 
However, validity of the perturbation theory implies upper bounds on the value of the 
charged Higgs-fermion Yukawa couplings. Requirements on the H + tb couplings yield an 
approximate region of phenomenological interest of 0.2 < tan/5 < 200 |3lJ. These limits 



16 




Figure 10: The search results, shown as 95% CL exclusion contours, in the [m H +, tan/3] 
plane for m t = 175 GeV and two it cross section calculations. 



are sensitive to m#+. 

We compute minimum and maximum (one-sided) limits on tan/?at 95% confidence 
level. Since tan (3 represents a physical parameter, it must have some discrete value 
and not a continuous distribution. It is therefore appropriate to provide minimum and 
maximum values for tan (3 separately. For a given set of [mt,a(ti)}, the maximum allowed 
values of tan/? can be plotted as a function of m H +, and, similarly, the minimum allowed 
values can be plotted. These form the boundaries of the regions in the [ra#+ , tan /3] plane 
that are excluded at 95% CL. 

The lightly shaded region in Fig. [lOj represents the region of [m#+ , tan (3] excluded 
for m t = 175 GeV and a (it) = 5.5 pb. The cross section of 5.5 pb is the one calculated 



by Berger et al, for m t = 175 GeV, using jj, = m t for the renormalization scale [14 



The reason for choosing this particular cross section is that it is the largest among the 
theoretical predictions at this mt, and therefore, for a given number of observed events, the 
most conservative in terms of setting exclusion limits on the existence of a charged Higgs. 
The darkly shaded region demarcates the additional excluded regions if one chooses the 
cross section calculated by Catani et al., {<j{ti) = 4.8 pb), for the same renormalization 
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For m t = 170 GeV we have a(ti) = 6.5 pb (Berger et al., \x = m t ). According to this 
cross section, given our lepton+jets signal and background acceptances, we are most likely 
to have seen 33 or 34 events if BR(tt — > W + bW~b) = 1. This means that n & s = 30 favors a 
non-zero contribution from channels in which at least one of the pair-produced tops decays 
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into a charged Higgs. Since the acceptance of the selection criteria for such channels is 
smaller than that for it — > W + bW~b, the larger the assumed it production cross section, 
the greater the need for a contribution from a charged Higgs to explain the observed 
deficit. A cross section smaller than that found in the standard analysis would yield a 
likelihood distribution with a single peak at tan (3 = ~ 6, where BR(tt — ► W + bW~b) 
has a maximum. As the assumed it production cross section is reduced, the likelihood 
distribution becomes narrower (and the area under the curve decreases), which results in 
a larger region of excluded on either side of the peak value. Like cr(ti), the value of m t 
also affects the derived range of excluded [m H +, tan j3}. 

It is important to note that these results are still preliminary. Two aspects of the 
analysis are being worked on and will probably change the results. The first one is taking 
into account another H + decay channel, H + — > W + bb, which becomes important at 
tan/?< 2 once one considers m H +> 130 GeV or so. Incorporating this additional mode is 
expected to somewhat weaken the lower limit on tan /3at large values of ra#+ . The second 
issue which needs further studies is the effect of increased widths of top and H + near the 
upper and lower boundaries of the range of tan /Sunder consideration. 



5 Future Prospects - Run 2 

The next Tevatron collider run (Run 2) is scheduled to begin in the year 2000 with 
increased luminosity and energy. D0 will increase its integrated luminosity compared to 
that accumulated during Run 1 by about a factor of 20 (to 2 fb _1 ). Combined with the 
increased Tevatron cm. energy (from 1.8 to 2.0 TeV) which will yield a 35% increase in 
the it production cross section and with detector improvements, the number of observed 
top quark events will increase by about a factor of 40. Details of the expected precision 
for top related measurements during Run 2, including some which were not done in Run 



1 due to the small statistics, are described in refs 35, 36 . 



6 Conclusions 

D0 has made direct measurements of the top quark mass from single lepton and from 
dilepton events using the entire data sample from the 1992-1996 running period with 
an integrated luminosity of 125 pb -1 . The top quark mass is determined to be m top = 
172.1 ± 7.1 GeV/c 2 . D0 has also measured the production cross section. Thirty 
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nine events were observed in eight different leptonic decay channels with an estimated 
background of 13.7 ± 2.2 events. In the all-jets mode 41 events were observed with an 
estimated background of 24.8±2.4 events. The tt production cross section is measured to 
be 5.9 ± 1.7 pb at m top = 172.1 GeV/c 2 . 
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